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Abstract. We consider the spatially homogeneous Boltzmann equation for inelastic hard- 
spheres (with constant restitution coefficient a £ (0, 1)) under the thermalization induced 
by a host medium with a fixed Maxwellian distribution. We prove uniqueness of the sta- 
tionary solution (with given mass) in the weakly inelastic regime; i.e., for any inelasticity 
parameter a 6 (a , 1), with some constructive ao 6 [0, 1). Our analysis is based on a 
perturbative argument which uses the knowledge of the stationary solution in the elastic 
limit and quantitative estimates of the convergence of stationary solutions as the inelastic- 
ity parameter goes to 1. In order to achieve this we give an accurate spectral analysis of the 
associated linearized collision operator in the elastic limit. Several qualitative properties 
of this unique steady state F a are also derived; in particular, we prove that F a is bounded 
from above and from below by two explicit universal (i.e. independent of a) Maxwellian 
distributions. 

Keywords: Boltzmann equation, inelastic hard spheres, granular gas, steady state, point- 
wise bounds, tail behavior. 
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1. Introduction 

1.1. Physical context: driven granular gases. Kinetic models for dilute granular flows 
are based, as well documented ||], on a Boltzmann equation in which collisions between 
hard-spheres particles are supposed to be inelastic, i.e. at each encounter a fraction of 
the kinetic energy is dissipated. Such a dissipation implies that, in absence of energy 
supply, inelastic hard spheres are cooling down and the energy continuously decreases in 
time. In particular, the corresponding dissipative Boltzmann equation admits only trivial 
equilibria. This is no longer the case if the spheres are forced to interact with an external 
agent (thermostat), in which case the energy supply may lead to a non-trivial steady state. 
For such driven system (in a space-homogeneous setting), the time evolution of the one- 
particle distribution function f(v,t), v E M. 3 , t > satisfies the following 

dtf = QUf,f) + G(f), (1.1) 

where Q a (f, f) is the inelastic quadratic Boltzmann collision operator (see next section for 
a precise definition), while G(f) models the forcing term. The parameter a e (0, 1) is the 
so-called "restitution coefficient", expressing the degree of inelasticity of binary collisions 
among grains and the purely elastic case is recovered when a = 1. 

There exist in the literature several possible physically meaningful choices for the forc- 
ing term Q in order to avoid the cooling of the granular gas. The most natural one is the 
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pure diffusion thermal bath for which the particles are subject to uncorrelated random 
accelerations between the collisions yielding to the diffusive operator 

A/, 

where \i > is a given parameter and A = A.„ the Laplacian in the velocity vari- 



able. For this model, introduced in [TO, the existence of a non-trivial equilibrium state 



has been obtained in [ |13[ ] while the uniqueness (in some weakly inelastic regime) and 



the linear/nonlinear stability of such a steady state has been proved in Q2l|]. Other 
fundamental examples of forcing terms are the thermal bath with linear friction [^]; 
02(f) = A Af + «;div(i)/) with several range of parameters k,X and where div is the 
divergence operator with respect to the velocity variable. A particular case of interest is 
the one related to the following anti-drift operator 

&(/) = -MdM«/), M>0. 



The existence of an equilibrium state for such a forcing term has been proved in [ pop 
and is related to the existence of self-similar solutions to the freely evolving Boltzmann 
equation. Such a steady state corresponds then to a self-similar profile (the so-called Ho- 
mogeneous Cooling State) and both its uniqueness (still in some weakly inelastic regime) 
and its stability have been derived in []20|], providing a rigorous proof to the Ernst-Brito 
conjecture [|T2|] for inelastic hard-spheres in the weakly inelastic regime. 



1.2. Description of the problem and main results. In this paper we address a problem 
similar to the aforementioned ones but with a forcing term of different nature. Namely 
we consider a situation in which the system of inelastic hard spheres is immersed into a 
thermal bath of particles so that the forcing term Q is given by a linear scattering operator 
describing inelastic collisions with the background medium. More explicitly, we shall 
assume in the present paper that the forcing operator Q is a linear Boltzmann collision 
operator of the form: 

Q(f) =: C(f) = Q e (f,M ) 
where Q e (-, ■) is a Boltzmann collision operator associated to the (fixed) restitution coef- 
ficient e G (0, 1] and Mo stands for the distribution function of the host fluid. We shall 
assume here that this host distribution is a given Maxwellian with unit mass, bulk veloc- 
ity uq and temperature ©o > 0: 

The precise definitions of both collision operators Q a (f,f) and £(/), with their weak 
forms and the relations between pre- and post-collision velocities, are given in Subsec- 
tion 2.1. 

The existence of smooth stationary solutions for the inelastic Boltzmann equation under 
the thermalization induced by a host-medium with a fixed distribution has been investi- 
gated by two of the authors, in collaboration with J. A. Carrillo, in we refer to this 
paper the references therein for more information about the physical relevance of such 
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thermal bath of particles. To be more precise, it has been proved in [Q] that, for any resti- 
tution coefficient a G (0, 1], there exists a non-trivial smooth stationary state F a ^ such 
that 

Qa(F a ,F a ) + C(F a ) = 0. (1.3) 

The proof of this existence result is based on a dynamic version of Tykhonov fixed point 
theorem and is achieved by controlling the L p -norms, the moments and the regularity of 
the solutions for the Cauchy problem Ql.ip. Moreover, using the analysis of the linear 
scattering operator C, for elastic nonlinear interactions (i.e. whenever a = 1) one can 
prove easily that there exists a unique solution with unit mass to the equation 

Q 1 (F,F)+C(F) = 0. (1.4) 

Moreover, this unique distribution is a Maxwellian M (v) with bulk velocity uq and explicit 
temperature G# < Go- The knowledge of the equilibrium solution in the elastic case 
a = 1 will be of paramount importance in our analysis of the steady state F a in the 
weakly inelastic regime a ~ 1. All these preliminary results are recalled in Section 2. 

Uniqueness and qualitative properties of the steady distribution are still open problems 
for a < 1, and these are the main subjects of the present paper. To be more precise, as far 
as uniqueness is concerned, we prove the following: 

Theorem 1.1. There exists a® G (0, 1] such that, for any g > 0, the set 

& a (g) = |f q G Ljj, F a ^ 0, F a solution to flPl ) with J F a (v) dv = gX 

reduces to a singleton where h\ is the set of integrable distributions with finite energy. In 
particular, for any a G (ojq, 1], such a steady state F a is radially symmetric and belongs to 



Several further qualitative properties of the steady state F a are also given in the paper. 
In particular, we are able to derive pointwise estimates for the steady state F a which are 
uniform with respect to the inelasticity parameter a: 

Theorem 1.2. There exist two Maxwellian distributions M_ and M (independent of a) such 
that 

M{v) < F a (v) ^ M(v) Vv G R 3 , Va G (0, 1). (1.5) 

For the upper bound, the strategy of proof is inspired by the comparison principle 
of [jl4p and uses some estimates of []lp . For the lower Maxwellian bound, the proof is much 
simpler than the ones yielding (non Maxwellian) pointwise lower bounds for the forcing 



terms Q\ and G3 (see Q20| |21D ) which rely on the spreading properties of the quadratic 
inelastic collision operator Q a . Our approach relies uniquely on the properties of the 
linear collision operator C and, more precisely, on the explicit integral representation of 
the gain part C + derived in |@J. 

More general theorems analyzing possible global stability properties of the stationary 
solution are planned as future work. 
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1.3. Strategy of proof and organization of the paper. Our strategy of proof is inspired 



by the strategies adopted in Refs. [pO, |2~T| for different kinds of forcing terms. However, 
the peculiarities of our linear scattering operator such as its lack of symmetry and the 
exchange of momentum between grains and background, will require in some points a 
completely different treatment with respect to previous works on analogous problems. 

1.3.1. Main difference with respect to other forcing terms. Let us spend a few words in 
explaining the key differences (which will also be emphasized throughout the paper): 

• The quadratic operator Q a (f,f) preserves mass and momentum and both the 



forcing terms Gi and G3 considered in Refs. p0| , |2lj ] also do so. Therefore, for 
both these forcing terms, the mass and momentum of a stationary solution can 
be prescribed. This is no more the case whenever the forcing term is the linear 
scattering £ which does not preserve momentum. 

Moreover, while the collisional operator Q a tends to cool down the gas — dissi- 
pating kinetic energy — the forcing terms Gi and G3 have the tendency to warm 
it up in some explicit way. Precisely, for any nonnegative distribution /, 

/ g x {f) \v\ 2 dv = 6fig f while / Q 3 (f) \v\ 2 dv = 2/i£ f 

where Qf = f R3 f(v) dv is the prescribed mass density of / and £f = f R3 \v\ 2 f(v) dv 
denotes its energy. It is unfortunately impossible to quantify the thermal contribu- 
tion of the linear scattering operator C in such a closed way: indeed, since we are 
dealing with a linear scattering operator associated to hard-spheres interactions, 
the thermal contribution f R3 £(/) \v | 2 dv involves moments of / up to third order. 
Finally, it is not possible in our case to use the fundamental scaling argument 



of [EQ, 21]. Precisely, for the forcing terms Q x and G 3 studied in [Eg, 21], seal 



ing arguments show that it is possible to choose fi > arbitrarily and this yields 



the authors of Q2(j, |2ip to choose fi = fi a so that, in the elastic limit a — > 1, the 
dissipation of kinetic energy will exactly be balanced by the forcing term. Such a 
scaling argument cannot be invoked for the linear scattering operator C and this 
is again related to the fact that we are dealing here with hard-spheres interac- 
tions. Notice that, if C were the linear Boltzmann operator associated to pseudo- 
Maxwellian molecules, the scaling argument and, more generally, the whole strat- 
egy of [ f20| , |2l| ] would apply almost directly (still assuming Q a to be associated to 
hard-spheres). 

1.3.2. General strategy. Let us now explain the main steps in our strategy of proof. It 
is essentially based on the knowledge of the elastic limit problem and on quantitative 
estimates of the difference between solutions to the original problem and the equilibrium 
state in the elastic limit. Introduce the linearized operator in the elastic limit (where we 
recall that, for a = 1, the unique steady state is an explicit Maxwellian M.) 

£> 1 (h) = Q 1 (M,h) + Q 1 (h,M)+£h (1.6) 
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Given a G (0, 1], let F a and G a belong to & a (g). One has 

Q a (F a , F a ) + C(F a ) = = Q a (G a , G a ) + C(G a ). 

Then 

Sfi(F a - G a ) = Qi(M,F a - G a ) + Q\(F a - G a ,M) - Q a (F a ,G a ) + Q a (G a , G a ). 
It is easy to recognize then that 



Sfi(F a - G a ) = [Qi{F a - G a ,M) - Q a (F a - G a ,M) 

+ \ Qi(M,F a - G a ) - Q a {M,F a - G a 

+ [Qa(F a - G a ,M - F a ) - Q a {M - G a ,F a - G 

Assume now that there exist two Banach spaces X and y (independent of a) such that 

\\Qi(h,M)-Q a (h,M)\\ x 

+ \\Qi(M,h)-Qa(M,h)\\x<v(<x)\\h\\y, Va€(0,l) (1.7) 
where lim Q ^i i](a) = and there exists C > such that 

\\Qa(h,g)\\ x + \\Q a (g,h)\\ x ^C\\g\\ y \\h\\y, VaG(0,l), (1.8) 

then, 

\\^ > 1 (F a -G a )\\ x < (r/(a) + C\\F a -M\\ y + 

+ C \\G a - M\\ y ) \\F a - G a \\ y Va€(0,l). 

If moreover there exists cq > such that 

\\&i(h)\\ x >co\\h\\y Vk [J 6 a (0)c^, (1.9) 

a6(0,l) 

then one sees that 

c \\F a - G a \\ y < 5(a) \\F a - G a \\ y Va G (0, 1) 

with 

5(a) = rj(a) + 2Cmax{||F Q - M\\ y , \\G a - M\\ y } . 
Therefore, if we are able to construct an explicit a® G (0, 1) such that 

F a ,G a G & a with a G (a , 1] => 5(a) < l/c (1.10) 

then 

F — G 

All the technical difficulty is then to determine X and y su ch th at (|L7t), ( |T8| ) and ( |T9| ) 



hold true and to prove that the 3^-norm is compatible with fll.lO[ ) 
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(1) The proof that ( ]T77p and ( ]1 .8p hold true will be straightforward on the basis of 



known estimates of the collision operator Q a . 
(2) Notice that ( |TT9t ) means that Jz^ : U ae (o,i) C A" is invertible and the 

proof of such a property relies on a careful spectral analysis of 



(3) Concerning now estimate (1.10), it consists in proving that 



lim sup \\F a — Ailly = 

More precisely it amounts to providing a quantitative estimate on the distance 
between F a and the Maxwellian M in the elastic limit a — > 1. This is the most 
technical part of the uniqueness result. 

To be able to complete the above program, one begins with deriving suitable a poste- 
riori estimates on the steady state that shall be useful in the sequel. In particular, after 
estimating the high-energy tails of the solution f(v, t) to the Boltzmann equation ( |1.1| ) 
uniformly with respect to the inelasticity parameter a, it is possible to prove that, for any 
a e (0, 1], the stationary solution F a admits an exponential tail of second-order. Moreover, 
we obtain uniform lower and upper bounds on the energy of F a and this yields a control 
of i7 fc -norms. 

To prove the points (2) and (3) of the above program, we derive the spectral properties 
of the linearized collision operator in the elastic limit ££\ given by fll.6Q . As already 
mentioned, this quantitative spectral analysis of Jzfi resorts to very recent results [ |l6| ] 
which allow to extend a spectral gap result from a smaller (typically Hilbert) space % to a 
larger (typically Banach) space X. We apply these recent abstract results to both the linear 
scattering operator C (whose spectral analysis in a weighted L 2 -space has been performed 
in [Q, [T8[ ]) and to the linearized operator 5£\, This will allow us to prove point (2) of the 



above program. Moreover, this spectral analysis will also allow to provide a quantitative 
estimate on the distance between F a and the Maxwellian M. in the quasi elastic limit 
a — > 1 (see point (3) above) . We wish to emphasize here the fact that, with our approach, 
we prove the convergence of F a to M as a — > 1 without knowing a priori that the energy 
£ p a converges to that £m of the Maxwellian M. This is a major difference with respect to 



the papers Q20| , |21[ ] where, for the reasons already explained, it was possible to write down 
a relatively simple equation (in closed form) satisfied by the difference £p Q — Em- This is 
not possible in the present situation since, again, £ is a scattering operator associated to 
hard-spheres interactions. 

1.3.3. Organization of the paper. After recalling the precise definitions of the Boltzmann 
operator Q a and the forcing term C, we give a precise simple proof of uniqueness of the 
equilibrium in the elastic case in Section 2. Then, Section 3 is devoted to the derivation 
of the a posteriori estimates on the steady state for general restitution coefficients. The 
uniform pointwise estimates (Theorem |1.2p are proved in Section 4 while Section 5 is 



devoted to the proof of Theorem |1.1| . In an Appendix of the paper, several estimates on 
C and ^fi are derived which turn out to be useful for the spectral analyis performed in 
Section 5. 
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2. Preliminary results 

2.1. The kinetic model. Given a constant restitution coefficient a £ (0, 1), one defines 
the bilinear Boltzmann operator Q a for inelastic interactions and hard-spheres by its ac- 
tion on test functions ip(v): 

Qa(f,g)(v) ip(v) dv = — [[[ f(v)g(w)\v-w\ U>{v') - ip(v)) dv dw da (2.1) 

47r JM.3 J R 3 J§2 

with v' = v + (\v — w\a — v + w). In particular, for any test function i/j = i^{v), one has 
the following weak form of the quadratic collision operator: 

Q a (f,f)(v)Hv)dv = \ [ [ f(v)f(w)\v-w\A a mv,w)dwdv, (2.2) 

A a [ip](y,w) = — / (tp(v f ) + 4>{w') - i>(v) - V>H) da 

4vr J§2 (2.3) 

= A+[ip](v, w) - A-[if>](v,w) 

and the post-collisional velocities (i/, «/) are given by 

1 + a 1 + a 
u = v H — (|g|(7 — g), w = w — — g), q = v — w. (2.4) 

In the same way for another constant restitution coefficient e G (0, 1), one defines the 
linear scattering operator C by its action on test functions: 

C{f)(v)^(v)dv= [ [ f(v)M (w)\v-w\JM(v,w)dwdv, (2.5) 



where 



J e [il>](v,w) = ±- I (#>*)-#>)) da = J+[iP](v,w)- J-[tP](v,w). (2.6) 



where 

4-7T 

with post-collisional velocities (v*, w*) 

1 + e 1 + e 
v * = v H — (|g|<7 — g), w* = w — (|g|cr — g), q = v-w. (2.7) 

For simplicity we shall assume in the paper that the total mass of the particles governed 
by / and that of .Mo are equal. Notice that 

C(f) = Qe(f,M ) 

and we shall adopt the convention that post (or pre-) collisional velocities associated to 
the coefficient a are denoted with prime symbol, while that associated to e are denoted 
with * symbol. We are interested in the stationary solution to the following Boltzmann 
equation: 

dtf(t,v) = Q a (f(t,');f(t,-))(v)+C(f)(t,v) t>0, f(0,v) = f (v). (2.8) 
We proved the following in 
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Theorem 2.1 (Existence of stationary solutions). For any restitution coefficient a e (0,1), 
there exists a nonnegative F a e L\ n LP, p e (1, oo) with unit mass and positive temperature 
such that 

Q a (F a ,F a ) + C(F a ) = 0. (2.9) 
Moreover, there exists a steady state which is radially symmetric and belongs to C c 



Remark 2.2. Notice that the existence of a radially symmetric stationary solution to ( |2.9P is 
not explicitly stated in where more general host distributions than Mo are considered. 
However, since the Maxwellian distribution Mq is radially symmetric, one easily checks that 
the property of being radially symmetric is stable along the flow of ( |2.8| ) (i.e. fo radially 
symmetric =^ f(t,v) radially symmetric for any t ^ 0) and therefore the fixed point 
argument used in [^] allows to build a radially symmetric steady solution to ( \2.9\ ). 

Notice that 

Q a (f, f) = Qt(f, f) ~ Q~(f, f) = Qt(f, f) ~ /£(/) 

where 

= (/*!• \)(v) = I fH\v - w\ dw. 

JTSL3 

Notice that £(/) does not depend on the restitution coefficient a G (0, 1]. In the same way, 

£(/)(«) = C + (f)(v) - £"(/)(«) = C + (f){v) - a{v)f{v) 

where 

a(v) = (M * | • \)(v) = / M (w)\v -w\dw. (2.10) 

2.2. A basic observation in the elastic case. We begin with a basic observation concern- 
ing the elastic case. Precisely, when the quadratic operator is that for elastic interactions, 
i.e. for a = 1, one can prove in a very direct way that the steady state solution to the above 
problem is unique. Precisely, the background forces the system to adopt a Maxwellian 
steady state (with density equal to 1): 

Theorem 2.3. The Maxwellian velocity distribution: 

"M=G^rM-^}' C2.ii) 



with 



B# = i^G (2.12) 



3-e 

is the unique solution with unit mass to the equation 

Qi(F,F)+C(F) = 0. (2.13) 
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Proof. It has been proved in [O] that C(M) = 0. Now, since M is a Maxwellian distri- 



bution, it is also well-known that Q\{M.,M) = and this proves that M is a solution 



to ( ]2.13[ ). To prove that it is the unique solution with unit mass, one proceeds in some 
formal way for the time being assuming that F decays sufficiently fast at infinity; we will 
see that it is actually the case in the following section. All the proof can then be made 



rigorous thanks to the subsequent Theorem |3.3| . For any distribution F(v) ^ solution to 

( 

\M{v) 



( |2.13p , let us multiply ( ]2.13p with log ( jMk ) anc ^ integrate with respect to v. One gets 



and it is well-known from [ ]T7| ] and [ |TT| ] that both the integrals in the above sum are 
nonpositive. Therefore, 



/ 

J'K 



Q^F,F){v)\og[^L\ dv = Q. 



M(v) 

Since M. is a Maxwellian distribution, it is a well-established fact that 

Qi (F, F) (v) log M (v) dv = 0. 
Consequently F is such that 

Q 1 (F,F)(v)logF(v)dv = 



and the classical Boltzmann H -Theorem [JTTJ] asserts that F is a given Maxwellian and 
Qi(F,F) = 0. Consequently, one has C(F) = and, from the uniqueness result [|l7], 



F = M. □ 

3. A POSTERIORI ESTIMATES 

3.1. High-energy tails for the steady solution. We are interested here in estimating 



the high-energy tails both of the solution f(t, v) to fl2.8p and of the stationary solutions 



to ( ]2.13[ ) through a weighted integral bound. Our approach is reminiscent to the work of 



[0] recently improved in a series of papers ||| Q [], [T4[ |J]. 

Definition 3.1. We say that the function f has an exponential tail of order s > if the 
following supremum 

r* =sup{r >0\T rs (f) := [ f(v) exp(r|u| s ) dv < +00} (3.1) 

JR3 

is positive and finite. 



We begin by showing that, for the solution to ( ]2.8| ), exponential tails of order s propa- 
gate with time if s E (0, 2]. The proof is adapted from several known results and follows 
the lines of d3l Section 6] . 
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Theorem 3.2. Let / be a nonnegative velocity function with f R3 fo(v) dv = 1. Assume that 
fo has an exponential tail of order s £ (0, 2], i.e. there exists tq > and s e (0, 2] such that 

fo(v) exp (r \v\ s ) dv < oo. 

Then, there exist < r ^ ro and C > (independent of a G (0, lp such that the solution 
f a {t, v) to the Boltzmann equation 

d t f(t,v) = Q a (f(t,-yj(t,-))(v)+C(f)(t,v) t>0, f(0,v) = f (v) (3.2) 

satisfies 

SU P / fa{t, v ) exp (r|t;| s ) dv ^ C < oo. (3.3) 

Proo/ We adapt the strategy of []8p following carefully the dynamical approach of [ p~4, [3p . 
For notations convenience, we shall drop the dependence on a for the solution to ( 3l2| ) 
and simply denote by /(t, its solution. Recall that, formally, 

/ f(t,v)exp(r\v\ s )dv=) —m S k(t) 

where 

m p (t)= f(t,v)\v\ 2p dv Vt^O.p^l. 



Therefore, to prove the result, it is sufficient to prove that there exists some r > (inde- 
pendent of a) such that 



OO I. 



^ T\ msk / 2 ^ conver g es for an y * > °- 

fc=0 

From the Cauchy-Hadamard formula giving the radius of convergence of a power series, 
it is enough to prove that, for any s € (0,2], there is some real number C = C(s) > 
(independent of t and of a) such that 

nisk{t)^k\C k Vt^O, VfceN. (3.4) 

2 

It is clear that, for any p ^ 1, the evolution of the p-moment m p (i) is given by 

d 



dt -mp(*) = Qp(*) + l pW 



where 



Q P (t)= [ Q a (f(t,-),f(t,-))(v)\v\ 2 Pdv and L p (t) = [ C(f)(t,v)\v\ 2p dv. 

Recall that the weak form of Q a and C are given in ( ]2.5p , ( |2.3[ ) and ( ]2.6[ ). Now, based 
upon a sharp version of Povzner's estimates, Bobylev, Gamba and Panferov [§, Lemma 1 
& Corollary 1] proved that, for anyp ^ 1, 

4[H 2!, ]K»)<7 a ,p(N 2 + HT 
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where, for any p > 1, 



7 



a,p 



1 fl + X 



h a (x) dx 



with h a (x) = I (g a (x) + g a (-x)) and 



9a{x) 



((1 - a)x + v / {l-a) 2 x 2 + Aa 



Vx E (-1,1). 



(1 + q)V(1 - «) 2 a; 2 + 4a 

Since we are looking for estimates which are uniform with respect to the inelasticity pa- 
rameter a, one notices that, as pointed out in JT9j 2^ 1, for any p ^ 1, 

4 



sup 7 Q) p < 7 P := min ( 1, 

«G(0,1) 



P+l 



In the same way, one obtains the following very rough estimate for J^~[\ ■ \ 2p ](v, w): 

J+[\ ■ \ 2p ](v, ™)^^J s2 (\v*\ 2p + K| 2p ) da ^ lp (\v\ 2 + \w\ 2 ) P . 
In particular, one obtains the following bounds: 

Aa[\ ■ \ 2P ](V, W) ^ 7P (\v\ 2 + H 2 Y ~ \V\ 2P - \ W \ 2P 

= -(1 - 7p ) (\v\ 2p + \w\ 2p ) + 7p ((M 2 + H 2 ) P ~ \v\ 2p ~ H 2p ) (3-5) 

and 

J e [\-\ 2p ](v,w) ^ lp (\ v \ 2 + \w\ 2 ) p -\v\ 2 *> 

= -{l- lp )\v\ 2p + J P [(\v\ 2 + \w\ 2 ) p - \v\ 2p - \w\ 2p ^j +j P \w\ 2p . (3.6) 
Then, as in [§, Lemma 2 & Eq. (4.5)], one sees that 

f + \w\ 2 y 



v — w 



\2V _ M 2p _ lwl 2 P 



fc=i 



^|2Cfe+l/2)|^|20—fc) + | w |2Cp-fc+l/2) |^|2^ (3 7:) 



where k r 



p+i 



using Eqs. (3.6)- 



is the integer part of Performing now the u and w integrations and 



( |37| ) we get 

/i 1) 



l pW < y £ (?) (^+i/ 2 (i)M p _ fc + m p _ fe+1/2 (t)M fc ) 
fc=i ^ ' 

+ 7 P (mi/2W^p + m (t)M p+1/2 ) 

f(t,v)Mo(w)\v - w\ \v\ 2p dvdw. 



l-7p 
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One estimates the loss term using Jensen's inequality (together with the fact that u = 0) 
to get 



/ 

Jm 



f(t,v)M (w)\v-w\\v\ 2p dvdw> / f(t,v)\v\ 2p+1 dv = m p+1/2 (t) 
and 

L P (t) < - -^m p+1/2 (t) + ^ (m 1/2 (t)M p + m (t)M p+1/2 ) + lp S p (t) (3.8) 

with 

~ 1 iv\ 

= 9 ( t ) ( m fc+i/2(0 M p-A,- + m p _ fc+1/2 (t)M fe ) . 
fc=l ^ ' 

In the same way, using Eqs. ( ]3.5p and ( |3.7p , the following estimate was derived in [|| 
Lemma 3]: 

Q p {t) ^ -(1 - lp )m p+l/2 {t) + lp S p {t) (3.9) 

where 

1 ^ / \ 

= 2 ( ) ( m fc+i/2 W«V^0) + rn p _ k+1/2 (t)m k (t)) . 
k=i ^ ' 

This yields to the following differential inequality for m p (t) : 

^ m p(*) < - 3(1 ~ 7p) m p+1/2 (t)+|^ (m V2 (t)M p + m (t)M p+1/2 )+7 P (^(i) + 

which is enough to prove that moments are uniformly propagated with time independently 
of a, i.e., for any p ^ 1, there exists C p > (independent of a) such that 

m p (0) < oo =^> supm p (i) ^ C p . 

Let us now introduce the renormalized moments 

, , rn n (t) 

*® := rWT6) ' Wlth a = 2/s ' 

where T(-) denotes the Gamma function and b > is a parameter to be fixed later on. 
Notice that, to get (^4), it suffices now to prove that, for any a ^ 1, one can find b > 
and some positive constant K > (both b and K independent of a) such that 

z p (t) ^K p Vp ^ 1. (3.10) 

An important simplification, first observed in consists in noticing that, for any a > 1 
and b > 0, 

S p (t) C T{ap + a/2 + 2b) Z p {t) 
where C = C(a,b) > does not depend on p and 

z p(t) = max {z k+1/2 (t) z p _ k (t),z k {t) z p _ k+1/2 (t)} . 
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In the same way, one proves easily that 

S p (t) ^ C T(ap + a/2 + 2b) Z p (t) for a ^ 1, b > 0, 

where 

Z p (t)= max {z k+1/2 (t) Cp-k, z p -k+i/2(t)(k} 

l^k^kp 

with (p = — ; ^ p . Then, using the approximation formula 



I>p + b) 



lim r ^ P + r) (a P y- r = l Va,i,r>0 
p->oo 1 (ap + e) 



together with the fact that j p = O \j\ as p — > oo, we get that, for sufficiently large p* ^ 1, 
there exist c\ , c 2 , C3 , C4 > such that 

^(t) + Cl p"/%(t) 1+1 / 2 P ^ C2 ^-i Cp + p -i+-/^ 1/a ) + C3p f +6-i^ (t) 

+ c 4 p a/2+fe ~ 1 £ p (i) Vr > 0, p ^ p*. (3.11) 

Now, since .Mo has an exponential tail of order 2, a fortiori it has an exponential tail of 
order s with < s < 2. Thus, for any a = 2/s ^ 1 and any b > 0, there exists C(a, 6) > 
and A > 1 such that 

M p ^C(a,b)T{ap + b)A p , Vp > 1, 
i.e. (p ^ C(a, 6)^4 P for any p ^ 1. Therefore, ( ]3.11[ ) becomes 

^(t) + d p a /%(t) 1+1 / 2p < C 2 (p" 1 ^ + p-i+-/2^+i/2^ + C3p i+&-i Zp (r) 

+ C 4 p2 i+6 ~ 1 2 p (t) Vt 5s 0, p > p* (3.12) 
for some positive constants Ci, C2, C3, C4 > where 

Z p (t) = max {z fc+ i/ 2 (t) ^~ fe , ^-fc+i/ 2 (i)^| • 

The key observation is that, for any p ^ p*, the functions Z p (t) and -H p (i) involve Zk(t) for 
^ p — 1/2 and do not involve z p (t). This is the reason why we will argue by induction in 
order to prove that, for any a ^ 1, if we choose < b < 1 it is possible to find K > large 
enough so that z p (t) ^ K p . First, because of the exponential integrability assumption 
on the initial datum /o, there exists Kq > such that z p (0) < Kq for any p ^ 1. Let us 
consider now po ^ p* > 1 such that 

2C 2 po 1 + (c 3 + a 4 )p - 1 < Ci, 

and let K > be such that 



X ^ < max sup Zk(t), Kq, 1, ^4 
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Since moments of f(t,v) are uniformly propagated, the existence of such a finite K is 
guaranteed. Defining now 

y p (t) ■= K p Vf ^ 

one can prove by induction (using also standard comparison of ODE's) that, for any p ^ po 
with 2p G N, y p (t) satisfies the differential inequality 

^E(i) + c x p a/2 y P (t) l+l ' 2p > C 2 (p~ 1 A p + p ^+^A p+l ' 2 ) 

+ c 3 p^ 1 z p (t) + c AP a i 2+b - 1 Z p (t), 

with moreover y p (0) ^ z p (0). One deduces from this that z p {t) < y P (t) = K p for any 
p ^ po and any t ^ 0. Notice that the comparison argument for ODE's is licit here since, 
again, for a given p ^ po, Z p (t) and Z p (t) involve only Zk(t) for k ^ p — 1/2. This yields 
the desired conclusion ( ]3.4p . □ 

We can now give a stationary version of the above Theorem in order to deduce the 
order of the exponential tail of the solutions to ( |2.13p . In the elastic case a = 1, as we 
already saw it, the solution to ( |2.13p is a Maxwellian and therefore has an exponential 
tail of order 2. For a given a & (0, 1), we look for the order of the exponential tail of 
the solution F a to ( |2.13p . Notice that, as shown in [Wj], the bounds obtained from Q a 
are actually uniform with respect to the coefficient a. This suggests that the order of the 
exponential tail of the solution F a shall be independent of a. Since, for a = 1, the order 
is s = 2, we infer that the solution F a to ( |2.13[ ) has an exponential tail of order 2. This is 
the object of the following Theorem: 

Theorem 3.3. There exist some constant A > and M > such that, for any a £ (0, 1] and 
any solution F a to ( |2.13p one has 

F a (v) exp (A\v\ 2 ) dv ^ M. 

Proof. The proof of Theorem p.3| follows exactly the same lines as that of Theorem |3T2 . 
We only sketch here the straightforward modifications. Recall that, for any r, s > and 
any a 6 (0,1], we defined 



*rAFa)= / 3 ^(E^H sfc )d, = ^^m f (a) (3.13) 

k=0 ' k=0 



where 



m p (a)= [ F a {v)\v\ 2p dv, p^O. (3.14) 
For any p ^ 0, we introduce now the following stationary moments 

Q P (a)= f Q a (F a ,F a )(v)\v\ 2p dv, L p (a) = f C(F a )\v\ 2p dv. 
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Of course, for any p ^ 0, Q p (a) + L p (a) = 0. Arguing exactly as above we get that 



3(1 - 7p )m; +1/2p («) < 7p (s p (a) + S p (a) + m 1/2 (a)M p + M p+1/2 



(3.15) 



where 



while 



S p( a ) = (fcj ( m fc+i/2(a) m p-fc(a) +m p _ fc+ i/ 2 (a)m fe (a)) 

^p( a ) = ( l) ( m fc+i/2(a)^ P -fc + m p-fe+i/2(a)^fc) • 
fc=i ^ ' 

To prove that the solution F a to fl2.13p has an exponential tail of order 2, as in the above 
proof (with s = 2) it is sufficient to prove that there exist C > and X > such that 

m p (a) < CT(p + l/2)X p , Vp > 1, Va G (0, 1]. (3.16) 

Notice that, since obviously TWo has an exponential tail of order 2, there exists Co > and 
Xq > 1 such that 

M p <C r(p+l/2)XP, Vp^l. 

Then, arguing as in the above proof, one gets that the above decrease of M p is enough to 
get (|3.16p by an induction argument as in the proof of Theorem [3.2|. □ 



Remark 3.4. The above Theorem provides some weighted L 1 space which contains all the 
stationary solutions for any a G (0, 1). Moreover, since the conclusion of the above result 
should hold for any a G (0,1), in particular, for a = 1 since F a = M\ is an explicit 

Maxwellian, one has A < A* with A^ := gQf • 

3.2. Uniform bound for the energy and control of the L 2 -norm. Upper bounds for 
the energy of the solution to ( J2.13 ) are easily obtained as a consequence of the above 
calculations or, more simply, from [|| Eq. (4.6)]: there exists E ni!a > such that 



\v\ 2 F n (v)dv < E n 



Va G (0, 1] 



where F a is a solution to fl2.13| ). In order to derive an uniform lower bound of E a (show- 
ing in particular that E a does not vanish in the elastic limit a — > 1), one shall actually 
derive an uniform lower bound of the L 2 -norm of any solution to (2.13): 



Theorem 3.5. Given a G (0, 1], any stationary solution F a to ( ]2.13p belongs to I? 
More precisely, there exists an uniform constant t 2 > such that 

As a consequence, there exists E m i n > such that 

E min ^ E a ^ E miiX VaG(0,l]. 



dv). 
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Proof. We prove the control of the I? norm as in [|20p. Precisely, let A > be fixed and 
let Aa(x) = "yX{x<a} + {Ax — ^f)X{x>A}> x The function is a C 1 -function over R 

2 

and lim^oo Aa(x) = \ for any x G R. In particular, for proving the claim, it is enough 
proving that there exists some positive constant c > not depending on a G (0, 1] such 
that 

limsup / A A (F a ) (v)dv < c. (3.17) 

A-^oo JR 3 

Let T A {x) := min (x, A) = A' A (x). Multiplying the identity ( |2.13p by T A (F a ) and integrat- 
ing over R 3 leads to 

f Q~(F a ,F a )T A (F a )dv+ [ C-{F a )T A (F a )dv 

JR 3 JR3 

= / Q+(F a ,F a )T A (F a )dv+ [ C + (F a )T A (F a )dv a€(0,l]. 

JR 3 JR3 

All the integrals in the above expression are nonnegative and in particular: 

I £-(F a )T A (F a )dv^ [ Qt(F a ,F a )T A (F a )dv+ [ C+(F a )T A (F a ) dv. 

JR 3 JR 3 Jr3 

Now, one estimates the left-hand side from below uniformly with respect to a as follows: 

f C-(F a )T A {F a )dv= [ F a (v)T A (F a )(v)dv [ M (w)\v-w\dw 

JR 3 JR 3 J R 3 

> c M [ F a {v)T A {F a )(v){\ + \v\)dv 
Jr 3 

where cm = inf„ f R3 M.q(w)\v — M dw/(l + \v\) is positive and finite. In particular, it does 
not depend on a. Then, as in p0[], since Aa(x) < xTa(x) one gets that 

CM 

Now, according to [ |20| , Step 2, Proposition 2.1], there exists 6 G (0, 1) such that, for any 
a G (0, 1], there is a constant C = C(E a ) > and A a > such that 

I Qt(F a ,F a )T A (F a )dv^C\\T A (F a )\\ 2 L ^- e) + [ A A (F a )(l + \v\) dv MA > A a . 

Jr 3 Z JR 3 

Notice that, though A a depends on the inelasticity parameter a, it will play no role since 
we are only considering the limit as A goes to infinity. Moreover, a careful reading of 
the proof of [ |20| , Prop. 2.1] shows that the constant C depends on a only through upper 
bounds of the energy E a . In particular, since we proved that E a ^ E m&x , one can set 
C = sup 0<a<1 C(E a ) < oo in the above inequality. One obtains finally 

[ A A (F a )(l + \v\)dv^C\\T A (F a )\\ 2 £- e) + [ C + (F a )T A (F a )dv MA > A a . 

* JR 3 JR 3 



I A A {F a ){l + \v\)dv^ I Qt(F a ,F a )T A {F a )dv+ [ C + (F a )T A (F a ) dv. 
Jr 3 Jr 3 Jr 3 
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One estimates now the last integral on the right-hand side owing to Theorem 1]. 
Precisely, according to the Cauchy-Schwarz inequality 

/ C + (F a )T A (F a )dv ^ \\C + (F a )\\ L 2\\T A (F a )\\ L 2 \/A> 

(notice that TA(F a ) £ L 2 for any fixed A since [1a(x)] 2 ^ Ax). Now, using the fact that 
£+(/) = Qt(fi-Mo), one deduces directly from [[j], Theorem 1] that 

\\C + (F a )\\ L2 ^C e \\F a \\ Ll \\M \\ Ll 

where C e > depends on the inelasticity parameter e. Since sup 0<Q<1 H-FaH^i < oo 
according to the result of the previous section, one gets that 

f C + (F a )T A (F a )dv^C \\T A (F a )\\ L 2 VA > 
Jr 3 

where Co > is a positive constant independent of a. Finally, we obtain 

^/ A A (F a )(l + \v\)dv^C\\T A (F a )\ftt e) +C \\T A (F a )\\ L 2 \/A > A a . 

* JR 3 

Since A A (F a ) ^ T A (F a ) 2 /2, this means that 

^\\T A {F a )\\l 2 ^ C\\T A (F a )fg- 0) + C \\T A (F a )\\ L2 VA > A a 

which clearly implies ( p.17 ). Now, it is a classical feature to deduce the uniform lower 
bound of the energy E a from the uniform control of H-FqII^ (see, e.g. Proposition 
4.6]). □ 



As in pQ , Prop. 2.1, step 8], a simple corollary of the above Theorem and the results 
of Section 6, Prop. 6.2] is the following uniform smoothness estimate: 

Corollary 3.6. For any k £ N, there exists Ck > such that 

\\ F a\\H k {R 3 ) ^ C k , Va G (0, 1]. 

In particular, there exists Coo > such that 

^Coo Va 6(0,1]. 



Proof. Notice that, from the uniform lower bound on the energy E a ^ E m \ n (a 6 (0, 1]), 
one notices that there exists some positive constant cq > such that 

E(F a )(v)= \v - w\F a (w) dw ^ c (1 + \v\) Vae(0,l]. 
Jr 3 

This, together with ||| Proposition 6.2], is enough to provide the uniform bound in 
H k (R 3 ) as in [ ]20| , Prop. 2.1, step 8]. Now, by Sobolev embedding theorem, one gets 
the second part of the corollary choosing simply k > 6. □ 
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4. Uniform pointwise estimates 

On the basis of the previous result, we derive in this section pointwise estimates for 
the steady state F a which are uniform with respect to the inelasticity parameter a. More 
precisely, we shall prove that there exists two Maxwellian distributions M and M (inde- 
pendent of a) such that 

M(v) < F a (v) < M(v) Vt> G R 3 . 
We will treat separately the upper bound and the lower bound. 

4.1. Uniform pointwise upper Maxwellian bound. The strategy of proof is inspired by 
the comparison principle of [[14}] and uses some estimates of [p. Precisely, a first general 
comparison principle is the following: 



Proposition 4.1. Let a £ (0, 1] and F a be a solution to ( 2.13| ) with unit mass. Assume there 



exists a measurable subset U o/]R 3 (with nonzero Lebesgue measure) and a measurable and 
nonnegative distribution G = G(v) such that 

Q a (G,F a ) + C(G) <0 for any v£U; (4.1) 

and 

F a (v)^G(v) for any weR 3 \«. (4.2) 
Then, F a (v) < G{v) for almost every v G R 3 . 



Proof. As already said, the proof follows the strategy of QT4] , Theorem 3] which is given 
for the time-dependent (space inhomogeneous) elastic Boltzmann equation. We adapt it 
in a simple way for granular gases in spatially homogeneous situations. One notices first 
that, for any nonnegative distribution g ^ and any distribution / 

/ Q a (f,g)(v)sign(f)(v)dv ^0. (4.3) 

Indeed, according to ( |2.1P , the above integral is equal to 

-f 

and the conclusion follows since <?(i>*) ^ while f(v) (sign(/) (i/) — sign(/(u))) ^ for 
any v, v* G M. 3 . For the same reason 

/ £(/)(u)sign(/(i;)) dv ^ for any distribution /. (4.4) 

Jm.3 



f(v)g(v*)\v - u*| (sign(/)(t; / ) - sign(/(u))) dvdv*da 



Now, after multiplying ( ]2.13p by sign(FQ, — G) and integrating over IR , one gets 

0= f (Q a (F a ,F a ) + £(F a )) sign(F a -G)dv= [ (Q a (G, F a ) + C(G)) sign(F Q — G) dv 

jR :i Jm 3 

+ f (Qa(F a - G, F a ) + C(F a - G)) sign(F a - G) dv 

JR3 
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and this last integral is nonpositive according to Q4.3P and ( |4.4| ). Therefore, 

/ (Q a {G,F a ) + C{G)) sign(F Q -G)d^0. 

Moreover, using the mass conservation property of both the collision operators, one can 
rewrite the above inequality as 

/ MG,F a )+C(Cf)) S!gnft ; G) + 1 dO0 
Jm.3 2 

where (sign(F a — G) + l)/2 is always nonnegative. We split this integral over U and its 
complementary. Whenever v £ U, by assumption ( |4.2p one has (sign(F Q (f;) — G(v)) + 
l)/2 = 0. Thus, the integral over M 3 reduces to the integral over U, i.e. 

/ (Q a (G,F a) + C(G)) Signft ; G) + 1 dO0. 

Now, according to our assumption ( |4.1p , the above is the integral of a nonpositive mea- 
surable distribution. Therefore, 

(Qa(G, F a ) + C{G)) (sign(F Q - G) + 1) = almost everywhere over U. 

Using again ( ]4.1p we get that sign(F a (t>) — G(v)) = — 1 for almost every v £ U which 
proves that F a (v) ^ G(v) for almost every v € R 3 . □ 

Now, in order to prove that every steady state F a is bounded from above by an universal 
Maxwellian distribution, we only have to determine a Maxwellian distribution G and a 
measurable subset U for which the above ( ]4.1| ) and ( ]4.2p hold true. We will need the 
following general result, proven in [Q, Proposition 11] that we state here for hard-spheres 
interactions only: 

Theorem 4.2 (Alonso et aZ.[Q]). Let 1 < p,q,r < oo with 1/p + 1/q = 1 + 1/r. Then, for 
a > there is a positive constant C a > such that 

llGiC/.^M- 1 ]^..^ < C7 a ll/M- 1 )!^^ IkM- 1 ^,^, Vae(0,l] (4.5) 
where M a (v) = exp(— o|i;| 2 ), v G M 3 . 



Remark 4.3. Notice that in [p the constant appearing in Q4-5p is actually given by C a Ci a , 
where C\ A is given by [[l], Eq. (6.10)] and depends only on a, while C a is given by [Q, 
Eq. (4.4)] and depends on the inelasticity parameter only through (1 — a) 2 . Bounding this 
last quantity simply by 1, one sees that sup a6 ( 01 ] C a < oo, thus obtaining a constant C a 
in ( |4.5p which does not depend on the inelasticity parameter a. 

This leads to the following 

Theorem 4.4. For any positive number a < min( 2,^,4.) where A > is given in Theo- 
rem [3.3| , there exists a positive constant p a > (independent of the inelasticity parameter a) 
such that 

F a (v) < exp ( - a\v\ 2 + p, a ) Vt;G]R 3 ,VaG (0,1]. 
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Proof. Let us fix a < mm(-^,A) and set M a (v) = exp(-a|v| 2 ), v G M 3 . As in [Q], one 
shall apply Proposition |4.1| with 

U = {v G M 3 , \v\ > R} 

for R > sufficiently large and with G(v) = K a M a (v) and K a to be determined. The 
technical part is to prove that Q4.1.P holds true for R > large enough. First, one has 

Q-(M a ,F a )(v) +£-(M )(u) = M a (v) ( [ F a (w)\v - w\ dw + / M (w)\v - w\ dw 



Recall that according to Theorem |3.3| , sup a6 ( 01 j f R3 F a (v)\v\ dv = m\ < oo. Therefore, 
since \v — w\ > \v\ — \w\ and both F a and Mq have unit mass, one has in a direct way 



Q-(M a ,F a )(v)+C-(M a )(v) > 2M a (v) \v 



mi + mo 



g 



(4.6) 



where m = / R3 Mq(w)\w\ dw. Now, to estimate C + (M a ) = Q+(M a , Mo), one applies 
( |4.5p with / = M a , g = Mo and (p,q,r) = (oo,l,oo). Since a < one sees that 
|| A^oM™ 1 Ili^RS) < 00 while trivially || MqM" 1 H^oo^a) = 1. Thus, there exists a positive 
constant c\ (a) > such that 

C + (M a )(v) SC ci(a)M a (t;) Vtj G M 3 . 



To estimate Q„(M a ,F Q ), one applies now ( \4.5\ ) with / = M a , g = F a and (p,q,r) = 
(oo,l,oo). Since sup a J" R3 F a (v) exp(^4|v| 2 ) dv < oo according to Theorem |3.3| , one sees 
that, for any a < A, 

SUp ll-FaM-^ljlma) < OO 
a€(0,l] 

and therefore there is a positive constant 02(a) > such that Q+(M a , F a )(u) ^ C2(a)M a (7;) 
V« £ I 3 . Gathering these two estimates, one gets the existence of a positive constant C a 
(independent of a) such that 

S+(M a , F a )(v) + C+(M a )(v) < C a M a (v) Vv G M 3 . (4.7) 

Combining ( |4~6[ ) and (]4~7|), one sees that, choosing R > " Tl+T ^ 0+C ' a , we have 

Qa(M«,Fa)(«) +£-(M a )(«) > C a M» ^ Q+(M a ,F a )(i;) + £+(M a )(¥) V|u| > R 

i.e. 

Qa(M ,F Q )(u)+£(M a )(u) ^0 WgW. 

Now, since there exists C > such that F a («) < C for any v G K 3 and any a G 
(0, 1] according to Corollary |3.6| , it is clear that one can find a positive constant K a = 
C exp(-ai? 2 ) > such that 

F a {v) SC if a M a (7;) \/\v\ < i?. 



With this choice of R and K a , the function G = K a M a (v) satisfies Q43f) and ( ]^2\ ) of 
Proposition |4.1| and we get our conclusion with fx a = log K a . □ 
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4.2. Uniform pointwise lower Maxwellian bound. We prove now a Maxwellian point- 
wise lower bound for the stationary solution F a which is uniform with respect to the 
inelasticity parameter. It turns out that the proof of such a result is much simpler than the 
ones yielding (non Maxwellian) pointwise lower bounds in the diffusively driven case [ ^Ij ] 
or for the homogeneous cooling state in J20[]. These two results rely on the spreading 
properties of the nonlinear inelastic collision operator Q a (in the spirit of similar results 



obtained in the elastic case in Q22Q). On the contrary, our approach relies uniquely on 
the properties of the linear collision operator £ and, more precisely, on the explicit inte- 
gral representation of C + derived in [ffl. We first prove a general lower bound for the 



time-dependent problem (2.8): 



Theorem 4.5. Let fo G L\be a nonnegative initial datum with unit mass and let f(t, v) be 



the associated solution to ( |2.8[ ). Then, for any to > 0, there exists a positive constant ciq > 



(which depends only on C > and to) such that 

f(t,v) ^ a exp(-7iH 2 ) Vv G R 3 ,Vi ^ to 
where 7l = 3±to! and ^ = 2^. 

Proof. The solution f(t,v) to ( |2.8p satisfies 

d t f(t,v) + &(J(t)){v)+a(v))f(t,v) = Qt(f(t), f(t))(v)+C + (/(*,-)) (v). 

Moreover, because of the propagation of moments uniformly with respect to a, there is 
some M\ > 0, independent of a such that 

sup / f(t,v)\v\dv ^ Mi < oo (4.8) 

so that there exists C2 > such that 

E(/(f))(«) + a{v) ^ c 2 (l + \v\) Vv G R 3 , W ^ 0. 
Therefore, the solution f(t,v) satisfies the following inequality: 

d t f(t,v) + c 2 (l + \v\)f(t,v)^ C + (f)(t,v) Vt^0,v£R 3 . (4.9) 
Now, according to the positive part C + admits the following integral representation 

£ + (f)(t,v) = [ k(v,w)f(t,w)dw 

where 

k(v,w) = C \v - w]' 1 exp { -j3 [ (1 + fi)\v - w\ + —. ^y- I } (4.10) 




with fj, = 2jy| ^ 0, Po = and Co > is a positive constant (depending on e and 0o). 
Moreover, the microscopic detailed balance law holds true 

k(v,w)M(w) = k(w,v)M(v) Vv,w £R 3 
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where M{v) is the Maxwellian distribution defined in (|2.3[): 



M{y) 



1 \ 



3/2 



2tt 



e# ) 



exp 



2Q* 



4(1 + M)A) 



Therefore £+(/)(*, u) =A^(u) / /c(w;,t;)/(t,w;)7W" 1 (w;)dw;. Since 

(H 2 -H 2 ) 2 



|v - wl 2 < 2|u| 2 + 2|w| 2 



straightforward computations yield 

Co 



and 



< 2\v\ 2 + 2\w\ 2 



v — w 



k(w, v) ^ 



\v\ + \w\ 



exp(— 7o|v| 2 ) exp(— 7i|w| 2 ) 



with 7 = 2/3 (l + /x + /j 2 ) and 7l = 2/3 (3 + 3/x + fi 2 ). Notice that A* 1 - 71 

any R ^ 2M 

f(t, w) dw ^ 



owing to the mass condition and fl4.8Q , for any R ^ 2M\ one has 

1 

2 

and 

Co 



-70 . Now, 
(4.11) 



exp 



(-(A»+7o)H 2 )x 



(27r e#)3/2 

/ ex.p(-~/i\w\ 2 )M~ 1 (w)f(t,w)- 
Jb(o,r) 

C exp(-(^ + 7o)|t;| 2 )x 



v\ + \w\ 



dw 



eM(A 9 -Ji)\w\ 2 )f(t,w) 

B(0,R) \V\ + H 

Hence, there exists Cr = Cq exp(— jqR 2 ) > independent of t ^ such that 



C + (f)(t,v)> 



\v\ +R 

This, together with ( $~% and ( ft.llQ , yields 



exp(-(^ + 7o)|u| 



B(0,R) 



f(t,w)dw. 



d t f(t,v)+c 2 (l + \v\)f(t,v) > exp(-7i|^| 2 ) V« G M 3 

2(Ju| + i?) 



from which we deduce 



/(*,«) £ 



2c 2 (M + i?) 2 
This clearly leads to the desired result. 



exp(- 7l M 2 ) (l - e -c2(i+M)A + e - C2 (i+k|)t /o ^) vt ^ 



□ 
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Remark 4.6. Notice that the above proof does not require the energy of f(t, v) to be bounded 
from below and the various constants involved depend only on the uniform upper bound on 
the first order moment fl4-8p . In particular, the lower bound of the previous Theorem 
shows that there exists a\ > 0, independent of a G (0, 1] such that 

inf / f a (t,v)\v\ 2 dv ^ ai > 



for any solution f a (t, v) to 

A stationary version of the above result is now straightforward: 
Theorem 4.7. There exists some positive constant ciq > such that, for any a G (0, 1], 

F a (v) ^ 1 exp(-a M 2 ) G K 3 . 

Proof. The proof follows the same paths of the previous one and is omitted here. Notice 
that the constant ao > does not depend on a G (0, 1) because the bounds provided by 
Theorem |3.5| are uniform with respect to the inelasticity parameter. □ 



The above uniform lower bound together with the regularity estimates of Corollary [ 
has important consequences on the entropy production. Precisely for any a G (0, 1], define 
the entropy dissipation functional, for any nonnegative g: 



@H, a (g) = 7—1 \v - w\g(v)g(w) 



8tt 



g{v')g{w') g{v')g{w') 
g(v)g(w) g(v)g(w) 



where the post-collisional velocities (v',w') = (v' a ,w' a ) are defined in fl2.4p . Notice that, 
for any nonnegative g for which all the integrals make sense, one has (see [13, |20| ] for 
details) : 

1-a 2 f 

Qa(d, d)(v) log g(v) dv = —S>H,a{g) H n— / g(v)g(w)\v - w\ dvdw (4.12) 

CT ./ro3 v ro3 



Then, arguing exactly as in [ ]20[ Corollary 3.4], we get the following 

Proposition 4.8. There exist ko and go G N large enough such that, for any ai > 0, there is 
some constant C > such that, for any g satisfying 

llffllif fe onii ^ a i> 9(v) ^ a 2 exp(-a 3 |t;| 2 ) 

one has 

Remark 4.9. Notice that, if /o G Ll Q n H ko is an initial distribution with unit mass, then, 
according to Proposition 6.3], the associated solution f(t, v) to Q2.8D satisfies 

sup\\f(t)\\ H k 0nL i <ai 
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for some positive constant a\ > 0. Hence, one deduces from Theorem ^5] and the above 
Proposition that there exists some constant C > such that 

sup |^,«(/(t)) - ®H,x (f(t))\ ^ C(l -a) Va G (0, 1]. 



5. Uniqueness of the steady state 

We aim now to prove that, for a G (0, 1) large enough, the steady state F a is unique, 
precisely, we show there is ao £ (0, 1) such that, for any g > and any a G (ao, 1), the set 

6 a ( g ) = i^F a G L\, F a ^ 0, F a solution to ( ]2.13| ) with j F a dv = £ j (5.1) 

reduces to a singleton. We first recall that Theorem |2.3| proves that it is the case in the 
elastic case: &i(g) = {gM} for any g > 0. On the basis of this easy result, we adopt 
the strategy described in the Introduction to prove the uniqueness of the steady state 
whenever a < 1. We begin by recalling the fundamental estimates of [ pTJ ] ensuring ( |1.7j ) 
and (Ot). 



5.1. Estimates on the collision operator. We recall here some results established in [ pip 
determining the function space in which the collision operator Q a depends continuously 
on the restitution coefficient a G (0, 1]. Let 

X = L 1 (?n~ 1 ) = ^(R^m^^dv), y = L\{m- 1 ) = L 1 (E 3 , dv) 

where 

m(v) = exp (— a\v\ s ) , a > 0, s G (0, 1]. 
Then, from [ ^T| , Proposition 3.2]: 

Proposition 5.1 (Mischler & Mouhot). For any a, a' G (0, 1) and any f G W^'^m.^ 1 ) and 

any g G L^m" 1 ), there holds 

\\Qi(f,g) - Q + Af,g)\U < p(« - <*')\\j 

and 

\\Qi(9J)-QU9,f)h<p(a-a')\\f\\ wl , lim _ 1) 

where p(r) is an explicit polynomial function with lim r ^, + p(r) = 0. 

With this proposition, one sees that (1.7) holds true for X = L 1 (m _1 ). Moreover, 
arguing as in [[l], Proposition 11], one proves easily that Q1.8P holds true: 

Proposition 5.2. There exists C > such that, /or any a G (0, 1) 

+ \\QZ(9,h)\\x < C||hb Vh,3 G 
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Proof. The proof follows from the very simple observation that 

\\Qi(h,g)\\x ^ WQiim-'h^n-'g^uw) (5.2) 



together with the well-known boundedness of the bilinear operator : L\ 
L}(M 3 ) ->■ L 1 (IR 3 ) (see, e.g. [[[], Theorem 1]). To prove $5% , one first notices that, 
for any h,g E X, one has 



\\Qi(h,g)\\x = \\Qi(h,g)m- 1 \\ L i= sup / Q+(fc, g)(v) (m"V) (v) dv. 

|H| ioo(K 3)=WR 3 

To estimate this last integral, one can assume without loss of generality that h,g,ip are 
non-negative. Then, using the weak formulation of Q+: 

Qa(h, g)(v ) (m -1 ^) (v) dv = / h(v)g(w)\v — w\ (m~ ip) (v*) dv dw 



where the post-collision velocity v* is defined by ( ]2.7p . Now, because of the dissipation of 
kinetic energy, since s e (0, 1], one has 

\v*\ s < (K| 2 + \w*\ 2 ) S ' 2 ^ {\v\ 2 + \w\ 2 ) s ' 2 < \v\ s + |^| s , 
i.e. m~ l (v*) < m~ 1 {v)m~ 1 {w). Therefore, 

Qa(h, g)(v ) [m~ l ijj) (v) dv ^ / (m~ ti) (v) (mT l g) (w)\v — w\ip(v*) dv dw. 



One recognizes that this last integral is equal to / Q^(m 1 h,m 1 g)(v)ip(v) dv and this 
proves (gg). □ 

5.2. Spectral properties of £ and Jzfi in A\ The spectral properties of both the linear 
Boltzmann operator C and the linearized operator «£? in % = L 2 (7W~ 1 ) are recalled in 
the appendix. In particular, it is well known that, for both these operators, is a simple 
eigenvalue associated to the eigenfunction M and both operators admit a positive spectral 
gap in %. We shall show that the same is true in the larger space X. To do so, we adopt 
the general strategy explained in the recent paper [ |16[ ] . First, we notice that H is a dense 
subspace of X. Moreover, if we denote, as in the appendix, Jz?2 as the linearized Boltzmann 
operator Qi(-,M) + Qi(M, •) + C in H, one has 

with jgfx : 0(jSfi) C * -> AT by 

JS?i(fc) = Qi(M,h) + Q 1 (h,M) + Ch, Vhe@(& 1 ) = y. 

In the same way, with the notations of the appendix, Cw = L- For the linear Boltzmann 
operator C, we have the following 
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Theorem 5.3. The spectrum of £ in X coincides with that of L in H. As a consequence, 
= span(A^) and £ admits a positive spectral gap v > 0. In particular, if 

x = {fex- f /cb = o}, y = {f € y ; [ fdv = o} 

then jY{£) n X = {0} and £ is invertible from y to X. 

Proof. As already mentioned, we adopt the general strategy explained in the recent pa- 
per [jl6p. Precisely one proves that £ splits as 

£ = A + B 

where 

(i) A : X ->■ U is bounded; 

(ii) the operator B : S'(B) —?■ X (with Si(B) = y) is a-dissipative for some positive 
a > 0, i.e. 

I s\gnf(v)Bf(v)m- l (v)dv^-a\\f\\ x V/ € y. (5.3) 
To do so, we use the estimates on £ derived in the appendix. For any R > 0, set 



Af(v) = £ + (x {H< K}/) («) = / «>)/M dw. 

Using Minkowski's integral inequality (with measures M^ 1 ^) dv and \ f(w)\ dw) one gets 
easily 

1 /2 

\\Af\\ n < [ \f(w)\( [ k 2 (v,w)M- 1 (v)dv) dw. 



'\w\^R 

Now, still with the notations of the appendix, one has 



k 2 (v,w)M' 1 (v)dv = M^iw) / G 2 (v,w)dv 



and, using Lemma A.|l|, with p = 2 and q = 0, there is some positive constant c 2 > such 
that 

f k 2 {v,w)M~ 1 {v)dv ^c 2 M~ l {w)(l + \w\y 1 ^^M^iw) Vw G I 3 . 

Thus, 

\\Af\\u<c[ \f{ W )\M- l ' 2 (w)dw 

J\w\^R 

and, since the domain of integration is bounded, there is some positive constant cr > 
such that || Af ^ 1 1 / 1 1 a* which proves point (i). Now, we prove that R > can be 
chosen in such a way that 

Bf{y) = £f(v) - Af(v) = £ + (x { \.\>R}f)(v) - a(v)f(v) 
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satisfies the above point (ii). For any / e y, set /(/) = f R3 sign f(v)Bf(v)m 1 (v) dv. One 



has 



r (f)= / sign/(v)£ + (x{|.|> J R}/)/(w)m 1 (v)dv- / a(v)\f(v)\m 1 (v)dv 



/ sign/(w)m 1 (u)du / k(v,w) dw — / cr(i;)|/(«)|m 1 (u)dv 

il 3 J{|w|>R} JR 3 

/ \f(w)\H(w)dw-a [ (l + luDI/CuJIm- 1 ^)^ 

J{|to|>/?} JR 3 



where we used the fact that cr(V) ^ <ro(l + H) for some positive constant oo > and set, 
as in the appendix, 



H{w)= / k(v,w)m 1 (v)dv, Vu> G M 3 . 
Jr 3 

Then, using Proposition A.|^, there is some positive constant .K" > such that 



1(f) \f(w)\(l + \w\ 1 - s )m- 1 (w)dw - a (1 + \v\)\ f (v)]™' 1 (v) dv . 

J{\w\>b} Jm.3 

In other words, 

1(f) < -ao ( \f(v)\m-\v) dv 



'{NCR} 

+ 



|/(u)| (K(l + l^l 1 " 8 ) - <r (l + M)) m~ l (v) dv. 

{\v\>R} 

We choose now R > such that if (1 + H 1_s ) — <to(1 + M) ^ — f° r all M > R (which 
can be done since s > 0), so that 



1(f) ^-a \f(v)\m(v)- 1 dv = -a \\f\\x, (5.4) 
Jr 3 

i.e. B satisfies ( ]5.3p with a = co- We conclude now with []16]. □ 

An analogous result holds true for the linearized Boltzmann operator S£ . Precisely, 

Theorem 5.4. The spectrum of ££\ in X coincides with that of Jzf 2 in % with Jf(££{) = 
span(M) and S£\ admits a positive spectral gap v > 0. In particular, jV(J£i) n X = {0} and 
jSfi is invertible from y to X. 

Proof. The proof follows exactly the same lines as the above Theorem ET3]. Precisely, one 
proves that Jzfi splits as Jz?i = stf + 88 where srf and ^ are such that srf : X — > H is 
bounded while the operator 83 : &(88) -> A? (with S>(88) = y) satisfies 

signf(v)mf(v)mT l (v) dv < -a||/|U V/ € y (5.5) 
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for some a > 0. One recalls that = Sff + Cf where Sff = Qi{f,M) + Qi(M, /). 
Now, it is well known that 



S?f = Sf+f- tTl (v)f(v)-M(v)[ / \v-w\f(w)dw 
where o\{v) = L a \v — w\M.(w) dw ^ a±(l + \v\) for some positive constant a\ > and 



jgf+/t = Q+(/i,7W) + Q+(A^,/i) = 2Qf(h,M). Then, it is easy to recognize (see e.g. 
that 



£>+h(v)= [ K 1 (v,w)h(w)dw 




with 

= Ci|t> — exp /3i ^|t> — w\ + 

where Ci > and (3\ = . In other words, J? + has exactly the same form of C + (with /3o 
replaced by /3i). In particular, Proposition A. ||| still holds if w) is replaced by _fTi(?;, w). 
One defines then 

= Sf+( X{l . l<R} f)(v) + £ + ( X{l . im fKv) 
(k(v, w) + K\(v , w)) f(w) dw. 

and s^2f(v) = —M(v) (J" R3 \v — w\f(w) dw) . It is clear that ^ is bounded from X to % 
with the very rough estimate H^/H^ ^ \\f\\ 2 x Jr3(1 + M) 2 -M( v ) dv . Moreover, with the 
same estimates as above (using the fact that the expression of Ki(v, w) is very similar to 
that of k(v, w)~), one proves that, for any R > 0, srf\ : X — > H is bounded. We define then 
£/ = srf\ + gfo so that 80 ■ X — > H is a bounded operator. Now, set 

&f{v) = Sfifiv) - *tf{v) = ^ + {x {V \> R} f) (v) + £ + (X{\.\>R}f)(v) - u(v)f(v) 

where v(v) = a{v) + ai(v) ^ i/*(l + \v\) with v* = cjq + a\. The estimates in the proof 



of the above Theorem |5]3] show then that there exists R > large enough such that 



satisfies (5.5) with a = z/*. We conclude as in [16]. □ 



5.3. The elastic limit a — )■ 1. The main result of this section provides a following quan- 
titative estimate on the distance between F a and the Maxwellian M. We recall the def- 



inition of & a {(>) in ( ]5.lD and, for simplicity, & a (l) shall be denoted <5 a . One has the 
following: 

Theorem 5.5. There exists an explicit function 771(a) such that lim Q ^i 771(a) = and such 
that for any cto G (0, 1] 

sup \\F a — A4\\y ^ 771(a) Vae(ao,l]- 
F a ee a 
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Remark 5.6. The fact that the above conclusion does not necessarily hold for inelasticity 
parameter a ~ is related to the estimate Q5.8P hereafter. Notice that this is no major 
restriction since the above result has to be interpreted as a result of uniform convergence to 
M whenever the inelasticity parameter a goes to 1 . 



Let us now come to the proof of Theorem ET5] which follows the paths of the corre 



sponding result in Q20]. Let M a denote the Maxwellian with the same mass, momentum 
and temperature as F a : 



M a (v) 



1 



27r6 a 



3/2 



exp 



v - u n 



2G a 



where 



u a = vF a (v) dv £ K 3 and Q a = - \ v - u a \ 2 F a (v) dv > 0. (5.6) 

One can prove the following result 

Proposition 5.7. Let a £ (0, 1] be fixed. For any q > and any 5 > 0, there is Cs(q) > 
such that the estimate 



\\F a - M a f L t S ^ C 5 (q)(l - a) VF Q £ & a , Va £ (a , 1]. 

9 

Proo/ Let a G (ao, 1] De fixed and let F Q e © a . The stationary solution F a satisfies 



(5.7) 



( ^ I 



(^1 



which, from [ ]T7| , Theorem 2.1], yields 

jf a Qa(F a ,F a ) log ( ^ > 0, Va G (0, 1] 



Now, 



/ Q Q (F Q ,F a )log(^)d«= / Q a (F a ,F a )logF a dv+-^ [ Q a {F a ,F a )\v\ 2 dv 

JR3 \MJ J R 3 2fc>» J R 3 



and, using (4.12) 



Qa^a, ^a) log F Q du = + 



a 2 Jrs 



F a (v)F a (w)\v — w\ dv dw 



while 



Q a (F a ,F a )\v\ 2 dv 



1-a* 



F a (v)F a (w)\v — w\ 3 dv dw. 
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Consequently, one has 



@H,a{ F <*) < 



1 - a' 



a- 



F a (v)F a (w)\v — w\ dv dw 



1 - a 2 
166» 



F a (v)F a (w)\v — w\ 3 dv dw 



1 - a 1 



i ? a (w)F Q (u;)|t>-u;|d?;dw; (5.8) 



From the estimate of the moments of F a , this last integral can be estimated from above 
by some positive constant K > independent of a € (0, 1]. In particular, for any fixed 
ao € (0, 1], there is Co > such that 

W«)^ c o(l-a) VoG(oo,l]- 

The above estimate, together with Proposition f£8[ implies the existence of some C\ > 
such that 



@H,i(F a )^Ci(l-a) VaG(ao,l]. 



(5.9) 



Recall that @h,i is the entropy dissipation functional associated to classical (elastic) in- 
teractions and has been studied intensively in [24]. In particular, using the estimates of 



Corollary |3.6| and Theorem |4.5| , one deduces from the op. cit. (see also Q2l| , Theorem 
3.5]) that, for any 5 > 0, there is Cg > such 



(5.10) 



Then, from ( |5.9P , we get that, for any 5 > 0, there exists Cg > such that 



|Fa-7Wa||ii < Ci(l-a)^ 



Va G (ao, 1). 



Now, using Theorems |3.3| and |3.5| , by a simple interpolation argument, we get the conclu- 
sion. □ 



An easy consequence of the above Theorem is the following where we recall that the 
space y has been defined in the previous section. 

Corollary 5.8. For any 5 > 0, there exists an explicit constant Cg > such that, for any 

a £ (0, 1] 

\\F a -M a \\y^C(l-a)*t5! \/F a e& a , Va€(a ,l]. 

Proof. The proof relies on a simple interpolation argument from Proposition ||]7| (with 
g = l) and Theorem |3.3| . Recall that y = L\{m~ l ) where m(v) = exp(— a\v\ s ) for some 
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fixed a > 0, s G (0, 2). For any a G (ao, 1] and any F a G 6 a , one has 

1/2 



l-Fa " Anally < l^a(w) - M a (v)\ (v) dv\ 



1/2 

|F a (v) - M a (v)\ (v) exp(2a|n| s ) dv 



Moreover, according to Theorem |3.3| (and since the energy E a of F a can be bounded from 
below and above independently of a G (0, 1)), there exist A > and M > such that 

/ \F a (v) — M a (v)\ exp(A\v\ 2 ) dv < M VqG(q ,1]. 

Since there exists c = c(a, s) such that (v) exp(2a|u| s ) ^ c(o, q, s) exp(A\v\ 2 ) for any v G M 3 , 
one gets the conclusion with C = y/M c(a,s)Cs(l) where Cg(l) is the constant appearing 



in Prop. |5/7[ □ 

With the above Corollary one gets the following: 
Lemma 5.9. There exist explicit constants C > and p > such that, for any ao G (0, 1], 

||/:(i ? a)IU = ||Q Q (i ? a ,i ? a )IU^C(l-a) p VF Q G6 a , Vae(ao,l]. 
Proo/ Let ao £ (0, 1] be fixed and let a G (ao, !]■ F° r an Y G &a, one has 

-£(F Q ) = Q a (F a ,F a ) = Q a (F a ,F a - M a ) + Q a (F a ,M a ) 

= Q a (F a ,F a - M a ) + Q a (F a - M a ,M a ) + Q a (M a ,M a ). 

Thus, 

ll^cOII* \\Q a (F a ,F a -M a )\\ x + \\Q a (F a -M a ,M a )\\ x + \\Q a (M a ,M a )\\x. 
Then, from Prop. |5.2|, there exists C > such that 



\\£(F a )\\ X < C ||F Q - M a ||y (H^b + \\M a \\y) + ||Qa(A<a, A*«)||*. 

Moreover, since M a is the Maxwellian with same first moments as F a , it is easy to see 
that || Anally depends only on the energy E a = J* R3 F a (v)\v | 2 dv. Thus, on the basis of the 
a posteriori estimates derived in Section 2, namely Theorem |3.5| , one gets easily that 

SUp (\\F a \\y + H-Mo-b) < OO. 
aS(0,l] 

Therefore, there exists a positive constant C2 > such that 

\\C(F a )\\ x ^ C 2 \\F a -M a \\y + \\Q a {M a ,Ma)\\x- 

Now, to estimate ||Q a (7W a , A / J Q )||,t', one only notices that, since M. a is a Maxwellian, one 
has Qi{M a ,M a ) = 0, i.e. 

\\Q a (M a ,M a )\\ x = \\Qa(M a ,M a ) - Qi(M a ,M a )\\x. 
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Therefore, one can apply Proposition |5j] to get the existence of some polynomial mapping 
r H> p(r) such that 

\\Qa(M a ,Ma)\\x < p(l ~ ") \\M a \\ w i,i (m -i } \\M a \\x 

where lim r ^o p( r ) = 0. Again, since the various norms of M a only depend on the en- 
ergy E a , we deduce from Theorem |3.5| that there exists a positive constant C3 such that 

\\Qa(M a ,M a )\\x < C 3 p(l -a) Vae (oo, 1]. 
Consequently, there exist two positive constants C2 , C3 > and some polynomial function 
r H> p(r) with lim r _>o p( r ) = such that 

< C 2 \\F a - M a \\y + C 3 p(l -a) Va G (a , !]■ 
We get the desired estimate using Corollary |5.8| . □ 



The above Lemma allows us to conclude the proof of Theorem 5.5 



Proof of Theorem p75\ For any a G (0, 1), set g a = —Q a (F a ,F a ), we get 

£{F a ) = C(F a -M)=g a 

with F a — A4 E y. Since C is invertible from y to X (with bounded inverse) according to 
Theorem |5.3| , there is some positive constant c > such that 

\\F a -M\\y = \\C~ l {g a )\\y < c\\g a \\ x . 

According to the above Lemma linia^i ||p Q ||^ = which yields the result. □ 

Remark 5.10. Notice that Theorem \5.5\ combined with Corollary |5.g| shows that 

lim \\M a -M\\y = 

a— 5>1 

with some explicit rate, where M a is the Maxwellian with same mass, momentum and tem- 
perature as F a . This implies in particular that 

lim u a = while lim Q a = 0* 



where u a and Q a are defined in ( |5.6[ ). Since C does not conserve momentum, it is not clear 
how to prove convergence of the first moments of F a towards those of M in a direct way. No- 
tice that, for the forcing terms considered in previous related works [ |20| , |2l[ ], the convergence 
of F a to M was, on the contrary, proved thanks to the convergence of the momentum and 
temperature. 

5.4. Uniqueness. With this in hands, as explained at the beginning of the section, one 
can state the following: 

Theorem 5.11. There exists a e (0, 1] such that, for any q > 0, the set 



& a ( e ) = < F a € L\, F a 0, F a solution to Q2.13p with / F a 



dv 



reduces to a singleton. In particular, for any a e (a.Q, 1], such a steady state F a is radially 
symmetric and belongs to C c 
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Proof. Our strategy to prove the uniqueness result has been explained in Section |5| and we 



have already shown that the estimates ( |1.7p , ( |1.8p and ( ]1.9[) h old true with our choice of 
X and y while condition fll.lO| ) holds thanks to Theorem |5.5| . □ 



Appendix: Properties of the linear and linearized Boltzmann operators 

A.l. Spectral analysis of the linearized operator in L 2 (A4 _1 ). We consider here the 
spectral analysis of the linearization of B(f, /) + £(/) around the Maxwellian state M. 
Precisely let H denote the Hilbert space L 2 (M~ l ) endowed with the inner space 

(f,9)n-= I f(v)g{v)M~ l (v)dv, Vf,geH 



and let S£i denote the following unbounded operator in %: 

& 2 {h) = L(h) + C(h), V/i G ®(Sf 2 ) 

where L(/i) = Q\(h,M) + Qi(M.,h) is the linearized operator of the classical Boltzmann 
operator Qi(-, •). The domain of ££ 2 in H is 

0(jSf 2 ) = L\(M~ l ) = {/ = f{y) ■ J \f{v)\ 2 M-\v) (1 + M 2 ) 1/2 dv < ooj . 

The spectral analysis of the linearized operator L in % is a well-known feature of the 
classical theory Boltzmann operator (see [ |TT , Chapter 7], [[jjl Chapter 3]). In particular, 
L is a nonnegative self-adjoint operator in H with 

(h,Lh) H < V7t € #(JSf 2 ) (A.1) 
and N(L) = span {M, v\ M, v 2 M.^v^M., \v \ 2 M } . Moreover, the spectral analysis of C in 



% has been performed in [Q] and made precise in Q18]]. Here again, £ is a nonnegative 
self-adjoint operator in 7i and there exists \i > such that 

- </,£(/)>* > Mll/-^ll! 2( M-i) V / g (A.2) 
with = J R 3 f(v)dv and where some quantitative estimates of the spectral gap /i are 



given in [ |18[ ] . In particular, iV(£) = span One deduces directly from ( |A..lp and ( |A.2| ) 
that 

-(f,Mf))n > Mil/ - e/^Hia^-ij V/ e #(JSf 2 ). 

In particular, 

AT(J&? 2 ) = span{7W}. 

Moreover, it is not difficult to resume the arguments of both [ ]l"5| ] and [Q] to prove that 
there exists some nonnegative measurable function v{y) such that 

where _£f 2 c is an integral operator, relatively compact with respect to the multiplication 
operator / h> vf. Therefore, the spectrum 6(^2) of Jz?2 is made of continuous (essential) 
spectrum {A e R; A ^ — ^0} where vq = mf veR 3 v{v) > and a decreasing sequence of 
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real eigenvalues with finite algebraic multiplicities whose unique possible cluster point is 
—uq. Moreover, the spectral gap ^ of Jz?2 



H2 :=min|A: —A G (-%0),-Ae 6(jSf 2 ) \ {0} 

satisfies the quantitative estimate fi2 ^ ^ where «. is the spectral gap of £ given in ( |A2p . 

A.2. Estimates on the linear operator C. We now establish several important estimates 
on the linear Boltzmann operator C. Precisely, we recall first the spectral properties of C 
in % where we recall that % = L 2 (M~ 1 dv). To distinguish the linear Boltzmann operator 
in % and in X, one shall denote by L the linear Boltzmann operator in %: L : &{h) C 
H -> H with 

9(L) = L\{M~ l ) = {/ = /(«) ; l/MI 2 ^- 1 ^) (1 + M 2 ) 1/2 dv < ooj 

and, 

Lf(v)= f k(v,w)f(w)dw-a(v)f(v), V/ G 0(L) 
where k(v, w) is given by Q4.10Q and cr(-) is defined in ( |2.10p and satisfies 

a(v)= k(v,w)dw^ a (l + \v\) Vt; G M 3 



with do > 0. Moreover, the spectral structure of L has been studied in [Q, [I8p and can be 
summarized in the following: 

Proposition A. 1. The spectrum <5(L) of the operator L inW is made of continuous (essen- 
tial) spectrum {A G K ; A < — v$\ where vq = inf„ eR 3 a(v) > and a decreasing sequence of 
real eigenvalues with finite algebraic multiplicities whose unique possible cluster point is —vq. 
Moreover, is an eigenvalue of L associated to M. and L admits a spectral gap > such 
that 

fi :=min{A : -A G (-z/ ,0),-A G 6(L) \ {0} j ^ 5il±£2 > Q 



4^5 

with 77 = a/2@o erf -1 (~) where erf" 1 denotes the inverse error function, erf _1 (i) ~ 0.4769. 

Notice that several properties of the kernel k(v, w) have been derived in [Q] in the spirit 
of [ ]T0| ] . Precisely, one has 

(w) = k(w, v)M («) V!),weR 3 xl 3 

and, setting 

G(v,w) = M" 1 / 2 ^)^^,^)^ 1 / 2 ^), v,w G M 3 x M 3 , 

one has G(u, 10) = G(w, v) and the following holds 

Lemma A. 1. For any < p < 3 and any g ^ 0, there exists C(p, q) > such that 

l G (»-«)l'7TTfn?< ri^il. - V "' 6KS ' 

(1 + (1 + 
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We shall exploit this estimate to derive the following more general one in which alge- 
braic weights are replaced by exponential weight. Namely, one proves the following: 

Proposition A. 2. Set m(v) = exp(— a\v \ s ), a > and s G (0, 1] and 

H(w) = / A;(?;,u;)m" 1 (u)d?;, w G R 3 . 



v — w 



Then, there exists a positive constant K = K(e, a, s) > such that 

H(w) ^ K(l + H 1-5 )™,- 1 ^) Vro € K 3 . 

Proof. Recall that k(v, w) is given by ( |4.10[ ). Taking into account that ^zffi 
2 \ lZZ \ ' w we ma y rewr ite Q4-10D as 

k(v,w) = C \v - w| -1 exp |-/3 ^(2 + fi)\v - w\ + 2 ^ - ^ ■ v?j j. (A.1) 

Performing the change of variables u = v — w and using spherical coordinates (with g = \u\ 
and g\w\y = u ■ w) one gets easily 

H(w)=2ttC [ F(g,y)dgdy 
J A 

with A = [0, oo) x [-1, 1] and 

F(Q,y) = £»exp <j-/3 f (2 + /J,)g + 2\w\y] + a( g 2 + \w\ 2 + 2g\w\y 



Split A into the two regions of integration: 

A 1 = {(g,y) e A;3\w\y^ -2g} and A 2 = A\A l . 
Notice first that, since y ^ 1 and s G (0, 1] 

exp (a{g 2 + \w\ 2 + 2g\w\y) s / 2S j ^ exp (a(g + \w\) s ) ^ exp(a^) ex.p(a\w\ s ) V(g, y) G A 
Moreover, since (2 + + 2|ii>|y ^ (/U + 2/3) f? for any y) G we have 

/ -^(^,y)d^dy ^ exp(a\w\ s ) / d^ / pexp (-/3 (2/3 + /i) V + a£ s ) dy 
M Jo J-i (A.2) 

^ Ci exp(a|u>| s ) = Cim~ 1 (w) 
since the integral is convergent. 

Let us estimate now the integral over A% which is more intricate. For any (g, y) G A 2 , 
one notices first that 

g 2 + \w\ 2 + 2g\w\y < \w\ 2 - g 2 /3 and g < {3/2)\w\, 
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so that 



f(q,v) dy dg 



A 2 



/•(3/2)|w| . , , rl . , 

J eexp (a (\w\ 2 - q 2 /3) s/ 2 J dg J ^ exp (-/3 ((2 + n)g + 2\w\y) 2 ) dy. (A.3) 



To carry out the y-integral, perform the change of variables z = (2 + p)g + 2\w\y to get 



/l /* OO 

^ exp (-/? ((2 + fi) e + 2\w\y) 2 ^) dy < _/ ex P 



-/3o^ 2 ) dz = ^ 



for some explicit C2 > 0. Plugging this in ( |A..3[ ) we obtain 

C r(3/2)|to| 



/ F(Q,y) dy dg ^ y- 2 - gexp (a (\w\ 2 - g 2 /3) s ) dg. 

Ja 2 \ w \ Jo v 7 



Setting now x 



w 



g 2 /3, we obtain 



A 2 2\w\ 



M 2 



l\w\' 2 /4 

We observe now that, for any r > 0, 



ex.p(ax s ^ 2 ) dx ^ 



3C 2 
2 M 



M 2 



exp 



ax' s/2 )dx. 



(A.4) 



[ T eMax s/2 ) dx < i- r i--/2 f —x s l 2 ~ l eMax s / 2 ) dx 
Jo as Jo 2 



2 i- s/2 f d 
— r 1 s/z / — exp 

as Jo dx 



(a// 2 ) dx < — r x ~ s l 2 exp(ar s / 2 
as 



Using this in (A.4) for r = \w\ we get 



F(g,y)dgdy ^ — 
A 2 as 



1 s exp (a|u;| s ) . 



Putting together flA.2Q and ( ]A.5| ) we finally obtain the result. 



(A.5) 



□ 



Remark A. 1. Notice that, whenever s = 1, the above Proposition actually asserts that 
H(w) ^ Cm~ 1 (w) for any w & M. 3 . Moreover, for any f & X = L 1 (M. 3 , m' 1 (v) dv) , one has 

\\£ + f\\x^ [ \f(w)\H(w)dw 

where C + is the restriction of L to X. In other words, for s = 1, we get that C + : X — > X 
is a bounded operator. This is reminiscent from [|l], Theorem 12] where exponential moment 
estimates for Q£(f,g) (with non Maxwellian weights) are derived. Notice that, in [Q, The- 
orem 12], an assumption of strict inelasticity (corresponding here to e < 1) was required 
which is not needed in the above Proposition. 
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